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Abstract. We study the Dirichlet problem for the non-local diffusion equation ut = 
J{u{x + z,t) — u{x, i)} d/x(z), where ^ is a. function and "u = 1^9 on 9f2 x (0, 00)" has 
to be understood in a non-classical sense. We prove existence and uniqueness results 
of solutions in this setting. Moreover, we prove that our solutions coincide with those 
obtained through the standard "vanishing viscosity method" , but show that a boundary 
layer occurs: the solution does not take the boundary data in the classical sense on 
a phenomenon related to the non-local character of the equation. Finally, we show that 
in a bounded domain, some regularization may occur, contrary to what happens in the 
whole space. 
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1. Introduction 

In this paper we investigate a suitable notion of Dirichlet Problems in a bounded set 
n C for a class of non-local equations, starting from remarks and ideas developed in 
[3] to cover a more general framework. 

To begin with, let us consider our model equation posed in the whole space : 
(1.1) —(x,t)= {u{x + z,t) - u{x,t)} dfi{z) , 

where /i G Ai{W^) is a finite measure. Such non-local diffusion operators arise in various 
problems like diffusion of population (see [S] and the references therein for more on this 
topic), Levy processes [HEllH], fractional Laplace diffusion [7], and more generally, integro- 
differential pde's (see [I] and the references therein). Note that in the fractional Laplacian 
case, fi{z) = 1^1^^^'*' with < A < 2, a case in which /x is not a finite measure that we 
shall not consider here (we refer to |3j for a general study of the Dirichlet problem in such 
situations). 

Now, if the equation has a rather clear sense when posed in the entire space M^, the 
question of what is a Dirichlet problem in this context becomes less obvious. Indeed, due 
to the nonlocal character of the equation, we need to assign the values of u not only on 
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the topological boundary dQ, but on the complement of Q (see more precise statments 
below). 

Let us mention that a viscosity approach of such Dirichlet problems can be found in [2j 
and [6j, but it concerns mainly the case of singular [i.e., non-integrable) measures. We 
shall prefer here a similar approach to that which can be found in [S] where we considered 
an homogenous Dirichlet problem with fi, a continuous bounded symetric function. But 
we want to go further in two directions: first to consider non-homogeneous problems (to 
be suitably defined), and second to include functions /i that may be singular, but yet 
integrable. 

Aims and Organization of the Paper 

The aim of this paper is thus to give a rather complete view of the case when /x > is 
an integrable function in M". 

We shall first define a suitable notion of what is a Dirichlet problem for our equation in 
Section 2, and then we investigate solvability of the problem in Section 3. More precisely, 
we prove that for any continuous and bounded initial data uq and "boundary" data (f, 
there exists a unique solution u, continuous and bounded in x [0, oo). We also prove 
that comparison holds in this class of solutions. 

Section 4 is devoted to study the vanishing viscosity method which is another way to 
try to solve the problem : we add a viscous term — eAn to the equation in order to force 
the boundary data to be taken, and then try to see what happens in the limit as e — > 0. 
Well, it turns out that we get indeed in the limit the unique solution obtained in Section 
3, but this solution is strictly positive on dfl whereas is zero on the boundary for any 
e > 0. Hence, we are facing a boundary layer formation which proves that in general the 
boundary data will not be taken in the usual sense. 

In Section 5 we investigate some fine regularity properties of the solutions. It is known 
that in the case Q = M", no regularization occurs for such non-local equations, in the 
sense that the solution at time t > cannot have a better regularity than that of uq (this 
was proven in [5]). But in the case of bounded domains, it appears that however some 
kind of smoothing does occur, provided /i is not too much concentrated: the modulus 
of continuity of the solution decreases exponentially as time increases (see Section 5 for 
precise statments). 

We will focus on some other qualitative aspects in Section 6. In particular, we prove 
that the equation enjoys an positivity property similar to the heat equation (if uq > 
and not identically zero, then u{t) is automatically positive in all Q). But in the present 
situation, this effect is not purely diffusive as for the Laplace operator, it is due to the 
nonlocal character of the equation. Moreover we show that solutions are immediately 
positive also on the boundary dQ, which is consistent with the layer formation in the 
vanishing viscosity method. 
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Finally we collect some other qualitative remarks, examples and possible extensions in 
the last Section. 

This note benefited a lot from discussions with Guy Barles and Ch. Georgelin, whom I 
warmly thank. 

2. Definition of the Dirichlet Problem 

We shall assume throughout this paper that the measure has a density (still denoted 
by /i) which is a nonnegative, integrable fonction of total mass /i(M") = 1. In fact, if the 
density is integrable, we can always assume that the total mass is 1 by invoking scaling 
arguments: simply put f (x, t) = u(x, Xt) where A is the total mass of /x. The set fl C 
will be open, bounded, connected and regular (to simplify). 

Now, as we noticed above, to give sense to a Dirichlet problem one needs to prescribe the 
"boundary" value not only on the topological frontier dQ, but in fact in the complement 
of ft or to be more precise, on the set where the measure ^ will read some information. 
So let be the translated measure defined by t~^ii{z) = fj,{x + z), let supp (//) denote 
the support of the measure /i and let us define the following subsets of : 



O'^ = U supp (r- V) , and a'^O = O'' \ O , 



In other words, Cl^^ is what we call the domain of influence of fl. In short, if you want 
to know how things evolve in fl, you need to read information in the bigger set Cl^^ D Q, 
because the operator will look at it. In this direction, the natural notion of frontier dfl 
has to be replaced by the extended frontier d'^Q, which is of course not a ncgligeable set 
since it is regular {i.e., the closure of an open set in M.^). Then in this setting, we may 
construct a general theory for the following Dirichlet problem : 

du 



{V) 



/ {u{x + z,t) — u{x,t)} dijL{z) in Q X (0, oo) , 

u{x, t) = (p{x, t) in d^fl X [0, oo) , 

u{x, 0) = Uo{x) in Q . 



The preceeding formulation requires the function u to be defined in fl'^ x [0,oo), ex- 
tended by in order to give sense to the integral term. But as we shall see below, u 
will not be continuous when crossing the boundary dfl so that practically, we shall often 
prefer a different formulation, keeping the function u to be defined only within Q x [0, oo). 
Let us denote by K^{u) the following non-local operator : for all {x, t) & fl x (0, oo), 

K^(u)(x,t) — / u{x + z,t) diji{z) + / ip(x + z,t) diji(z) — u{x,t) . 



Thus we arrive at the following definition 
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Definition 2.1. Let fi be a finite measure, (f G C(9'^f2 x [0,oo)) and Uq G C{Q) with 
Uq = ip in dfl X {0}. By a solution of (V), we mean a function u G C(fi x [0, oo[) such 
that Ut = K^{u) in Q X (0, oo) in the sense of distributions and u{0) = Uq in Q. 

Note that writing ut = K^p{u) implicitly means that we are trying to solve the Dirichlet 
problem with "m = in d^Vt x (0, oo), while u is only defined in f2 x [0, oo). Note also 
that since we consider continuous solutions, Ut will be continuous and the equation will 
hold in the classical sense. 

A more general definition could be given, only imposing that the translations of (p and 
Mo be integrable with respect to yU. Thus, in the case when /i is a bounded function, it 
is enough to consider that ip and u are integrable with respect to the Lebesgue measure 
and the formulation makes sense. But in this paper, we shall restrict ourselves to the case 
/i G L^(M") (see next Section). 

Notice that if the support of /i is not compact, then we may have possibly Vt^ = M.^ 
and d'^Q = M^\i7, the case in which every point, even outside Q will infiuence the density 
at any other point x in Q. On the contrary, if /i is compactly supported, then the values 
of u outside Cl'^ will never enter into play, as soon as we look at the problem in Q. Of 
course the homogeneous case treated in ^ appears as a special case of this framework 
when ip = and fie C(M"). 

We may now come back to our first intent and make precise the formulation : if m is a 
solution in the sense of Definition 12.11 we can extend it as follows 



and for instance u = outside x (0, oo) (the values there are not involved in the 
integral). Then u satisfies : 



In this Section, we shall derive here a general existence and uniqueness result when n 
is a measure. In the following, the notation Cb{Y) means the set of continuous and 
bounded functions on Y, and let 



Theorem 3.1. Let /i G L^{R). Then for any Uq G Cf,{Q) and ip G Cb{d^Q x [0, oo)) with 
Uq = ip in dVt X {0}, there exists a unique solution u & X of the Dirichlet problem (V). 

Proof. This result is obtained by a standard Banach fixed-point theorem in X as follows: 
consider for to > fixed the Banach space Xt^ = C([0,to] x ^) and the operator T„g : 



(2.2) 





3. A General Theory of Existence and Uniqueness 



X = Cb{n X [0,oo)). 
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Xtg Xtg defined by 



/ K^{u){x, s) ds . 
Jo 



Let us show that Tug (u) remains in Xt^ , and notice first that boundedness is obvious since 
/i(R'^) = 1 implies that \\K^{u) {-,1)11^ < 2||m(-, t) ||oo. 

Now, take a;„ G ^ x G and t E [0,to]- Since u is continuous up to the 

boundary, to pass to the limit and check that Tu^iu) is continuous, we need only to 
investigate the term 

u{Xn + Z, S) d^{z)(\.S + I I (p{Xn + Z, S) (\^{z)(\.S . 



To treat both terms at the same time, consider the extension it defined in f l2.2p . Then 
u{xn + z,s) u{x + z, s) for any z such that {z + x) ^ dfl (recall that u and ip are 
continuous on either side of the boundary of Q), so that u{xn + ■)—>■ u{x + ■) almost 
everywhere in R" and u is uniformly bounded. By dominated convergence we get that 

u{xn + z)d^{z) / u{x + z)dfi{z) , 
Jo 

which proves tlicit ^?~u0 

(u) is continuous in all [0,to] x ^• 
Now it is easy to check that (recall again that yu(]R") = 1) 

sup \Tuo{v) - Tuo{w)\ < 2to ■ sup \v -w\, 

Qx[0,to] rix[0,to] 

which implies that is a strict contraction in X^^ provided to < 1/2. Hence in this case, 
we may use Banach's fixed point theorem in Xt^ to obtain a unique solution of u = Tuq{u) 
in Xjq. Once a solution is constructed on [0,^0], it remains to iterate the procedure and 
finally get a unique bounded, continuous solution in x [0, cxd[, which ends the proof. □ 

Proposition 3.1. For any solution u E X of {V), the equation Ut = K^{u) holds at any 
boundary point {x,t) G dfl x [0, 00), and also at t = in Q (with continuous extension of 
dtu at t = 0). 

Proof. This is done exactly as in the Theorem above : if (x„,t„) ^ (x,t) G x [0, cxd), 
then we can pass to the limit in the different terms of the integral equation f 1 1.11) so that 
we get the same equation at {x, t), even if it is a boundary point x G dQ, or t = 0. □ 

Proposition 3.2. Let ^ > be a L^iW^) function and assume that u E X is a solution 
of {V) with initial data Uq and boundary data (f ; v E X is a solution with initial data Vq 
and boundary data ip. Then if Uq < vq and (f < ip, we have u < v in all Cl x [0, oo[. 

Proof. Let Wc{x, t) = v{x, t) — u{x, t) + c/(T — t) defined in x [0, T], for some arbitrary 
c, T > 0. The continuous function Wc reaches its min at some point (xo,to) G x [0,T] 
and obviously to cannot be T. If to = then we are done, so it remains to investigate the 
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case when xq & Q and t > 0. Using the fact that necessarily the equation holds at such a 
point (even if it is located at the boundary), one gets the following : 

dtWc{xo,to) = = / {wc{xo + z,to) -Wc{xo,to)}dfi{z) + , 

Jr" [T-ty 

where Wc = v — u + c/{T ~t) (natural extensions of u and v defined by (12.21) ). Remember 
first that Wc = ip — ^ + c{T — t) > in d^Q x [0, T], and second that Wc reaches its min 
in X [0,T] at (xo,to)) so we arrive at 

c 



Wcixo,to) = / Wcixo + z,to) dfx^z) + 



(T - ty 

> I Wc{xo + z,to)dfi{z) > + 



'{{xo+z)en} 
> aWc{xo,to) + 



(T - ty 



(T - ty ' 

for some a G [0, 1[ (recall that we assume /i(M") = 1 and note that the case a = 1 leads 
directly to a contradiction since c > 0). Hence Wc{xo,to) > and passing to the limit as 
c — i> yields the result : u > v in Cl x [0,T]. Since T > is arbitrary, the comparison 
result holds. □ 



4. The Vanishing Viscosity Approach 



In this section we try another approach to construct a solution, namely using the 
vanishing viscosity method. It consists in adding a viscous term —eAu to the equation 
to get a classical solution that does take the boundary data, and then pass to the limit 
as £ — * 0. We shall do this in the case = for simplicity, the modifications when 
is continuous being rather straightforward. 

Recall that in the case (/? = 0, we denote by Kq{u) the non-local operator associated 
with the zero boundary data : 

Kq{u){x) = / u{x + z) d^{z) — u{x) . 

Proposition 4.1. Let e > 0. Then for any uq G Cq{Q), there exists a unique solution 
u = Us in C{Cl X [0, oo)) of the following problem : 

{ut — eAu = Ko{u) in X (0, oo) 
u{x,t) = on dnx{0,oo) 

u{x, 0) = Uq{x) in 

Proof. Existence of a solution follows from standard arguments once we have a comparison 
result and some bounds. In fact the comparison principle for this viscous equation is clear: 
it is obtained like in the case e = (see previous section) with standard arguments since 
the viscous term does not change anything here. Moreover, positive constants are solutions 
(with positive boundary values) so that for uq G Co(f2) we obtain a solution u which is 
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bounded by HmqIIoo- Thus u satisfies Ut — Au = f with / bounded, which imphes that the 
zero boundary data is taken in the classical sense, and that u G 0(0, x [0, oo)) fi C^(f2 x 
(0, oo)). Uniqueness follows also from the comparison principle. □ 

Theorem 4.1. As e ^ 0, converges pointwise and in L°°{Q x {0, oo)) -weak-k to the 
unique solution u G 0(0 x [0, oo)) of the following problem : 

Ut = Kq{u) in X (0, oo) , 
0) = Uo{x) in . 

Proof. Let us begin by using a weak formulation for the viscous problem : for any (j) G 
Cl{Q X (0, oo)) {i.e., regular and compactly supported in x (0, oo)), we have 

/ Us{t)(j){t) — / Uo0(O) — £ / / MeA0dxds= / / Ko{us)4>dxds . 
Jr" Jr" Jo Jr" Jo Jr"- 

Then since the sequence {ue)e is bounded in L°° as e — > 0, we can extract to get a sequence 
{us„) converging in L°°-wesik-k to some function u^,, bounded in fi x (0, oo). This allows 
us to pass to the limit in the weak fomulation: 

u^(j){t) - / Mo0(O) = / / i^o(M*)0da;ds . 







This equality proves that at least in the sense of distributions in fi x (0, oo), the bounded 
function u^. satisfies dtu^. = Kq{u^). Now, let us show that Kq{u^) is a continuous function 
in fi X (0, oo), so that the equation holds in the classical sense. Since m* is bounded and 
/i G L^(]R"), we can use the density of continuous functions in to prove that 



(x, t) + z, t) dfi{z) is continuous in fi x (0, oo) . 

J{(x+z)eQ} 

Indeed, this is just using the dominated convergence theorem for a continuous integrand 
as in the proof of Theorem 13. 11 and then passing to the limit in by a density argument. 
Then, since Ut is a continous function in fi x (0, oo), we can integrate in time the equation 
satisfied by v{x,t) = e^u{x,t) to get : 

v{x,t) — v{x,0) = / / v{x + z, s) dfi{z)ds , 

Jo J{(x+z)eQ} 

hence v is continuous and u^, also, in all fi x (0,oo). Finally, as we already did in the 
previous section, we can pass to the limit as x goes to the boundary to conclude that 
the equation is satisfied in fi x (0,oo). Moreover, u^, can be extended to get a solution 
u G C(fi X [0, oo)), so that u is the unique solution of the problem (by Theorem 13.11) . 

Let us now prove that the limit holds pointwise : take T > fixed. Since u^^ —>■ u in 
L°°-weak -k, we can always assume that the limit also holds in L^-weak. Let us make some 
energy estimates in for the solution : since u^^ = on the boundary, 

\Uef{t)+en [ [ \DUef= f \uo\' (t) + [ [ MUeJUe^ , 

n Jo Jn Jn Jo Jn 
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SO that in the hmit we get 

/ < / |moP+ / / Kq{u^)u^. 
Jn Jn Jo Jn 

But since this inequahty is in fact an equahty (obtained by multiphng the hmit equation 
by Ut, and integrating over Q x (0,T)), we get that over all Q x (0,T), 

lim ||-U£„||i2(f^x{o,T)) = ||'"*||L2(nx(o,T)) • 

So, since we have both the weak convergence and the convergence in the L^-norm, 
it turns out that the convergence holds in strong. Thus, up to another extraction, 
u^^ — M^, pointwise in x (0, oo) (T > is arbitrary). Finally, the arguments above show 
that for any other extraction of the sequence (we), the limit is necessarily u^, which implies 
that in fact it is all the sequence (us) that converges to m*, pointwise in f2 x (0, oo). Thus, 
Us ^ u in Q X [0, oo) since u is the continuous extension of to x [0, oo). □ 



5. Regularity Questions 

Let us first mention that obviously if ;U G L^(M"'), any continous solution u{x, t) is in 
time, which follows directly from the integral formulation: Ut is continuous in x [0, oo). 

Now, as far as the space regularity is concerned, we begin recalling that the equation 
does not have any regularizing effect when /i G L^. More precisely, it was shown in [5j in 
the case = M" (and continuous) that if u is the solution with initial data mq, then 
u{t) = e'^uo + v(t) where v is smooth, so that the regularity of uq is preserved for all 
times. See also [Ij for estimates in BUG in (in a more general context). 

We will prove that we can control the space modulus of continuity. If u G C{Q x [0, oo)), 
let us denote by 

uj{r],t)= sup \u{x,t) - u{y,t)\ , 

\x-y\<v 

and let us begin with a result in the entire space M": 

Theorem 5.1. Let = M", /i G L-^{W"-) and u G C(M" x [0, oo)), a solution with initial 
data uq having a modulus of continuity uq- Then we have : for any rj > and t > 0, 

uj{r],t) < uJo{r]) . 

Proof. Let us integrate in time the equation for v{x,t) = e^u{x,t) and denote by Uy the 
modulus of continuity for v : 

v{x,t) — v{y,t) = v{x,0) — v{y,0) + / / {f (x + 2;, s) — y(?/ + z, s)} d/i(z)ds , 

Jo Jr" 

so that if |x — ?/| < r], we can estimate the modulus of continuity as follows (recall that 
by assumption, /i(M") = 1) : 

\v{x,t) -v{y,t)\ <uJo{ri) + / u;^(?7, s)ds. 

Jo 
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Hence using Gronwall's lemma we get that u!y{r], t) < u!o{r))e^, and going back to u{x, t) = 
e~^v{x,t) gives the result. □ 

In this generality, the Theorem is optimal : take N = 1 and for instance iJ,{z) = 1/2 
on [—1; 1], fi{z) = if \z\ > 1. Then if ^0(2^) = ox, the solution is stationnary and so the 
modulus of continuity is constant. 

The situation is a little bit more complicated in a bounded domain Q since the boundary 
terms may induce some sort of discontinuity. So let us first restrict ourselves to the case 
(p — and let us introduce some notation : 

For any x e IR", let t'^D, ^ {z : {x + z) e Q} and for any two sets A, B,let AAB ^ 
{A \ B) \J {B \ A). Now, for any rj > 0, let us define the following quantities : 

X{7])= sup fi{T-^QnT-^Q) , 

(x,y)en2 

k-!/|<J) 

7(77)= sup //(r-ifiAr-^n). 

(x,y)en2 
k-!/l<') 

Notice that obviously, X{'r]):l{v) ^ [0; 1] ^^'^ that 7(77) goes to zero as 77 — > 0. Roughly 
speaking, X{r]) is related to the total mass of // in Q = (J t~^Q, while 7(77) is related to 

the mass of n involved at the boundary. 

Theorem 5.2. Let us assume that Q is bounded and regular, fi > is a measure, 
= and that u is a solution with initial data uq G Co(r2). Then the following estimates 
hold : 

• // X{r]) < 1, then for any t > 



uj{r],t)< (^a;o(?7) +7(?7) • Ikolloo^-^yp— je 



{X{v)-l)t 



• // X{r]) = 1, then for any t > 0, 

uj{r), t) < (uJo{r)) + -f{r)) ■ \\uq\\^ ■ . 

Proof. The proof proceeds with the same ideas as for the case Q = M" : integrating in time 
the equation for v{x,t) = e^u{x,t) and again denoting by cUy the modulus of continuity 
for V one gets : 

v{x,t) — v{y,t) — I I u(x + 2;, s) d//(2;)ds — / / ^(t/ + 2;, s) d//(2;)ds , 

SO that if I x — y I < 77, we can estimate the modulus of continuity as follows : 

\v{x,t) — v{y,t)\ < I I \v{x + z, s) — v{y + z, s)\(iiJi{z)(is + R{x,y,t) 

Jo Jr-'^QnT-'^n 
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where R{x,y,t) contains the integral over (r^^flATy^fl), that we estimate uniformly 
using the fact that 11^(5)1100 ^ ||'i*o||oo " '■ 

R{x,y,t) <\\uo\\oo-l{v) ■ eMs. 

Jo 

Thus, going back to the estimate, we have for any \x — y\ < rj, 

\v{x,t)-v{y,t)\<ujy{r),0)+ / {X{r)) ■ ujy{r], s) + -f{r)) ■ \\uo\\^ ■ e^} ds . 

Jo 

Passing to the supremum over all {x, y) such that — ?/| < 77, and using Gronwall's lemma 
we obtain: 

O0v{v,t) < (o;,(r/,0) +7(^) ■ ||%||oo ^ _ ^(^^^ 

and finally going back to u{x,t) = e^^v{x,t) yields the desired estimate. When A(?7) = 1, 
then the estimate is easier to obtain, using also Gronwall's lemma. □ 

Before stating a corollary, let us give a heuristic idea of the meaning of those estimates: 
if A = 1, we see that the estimate on the modulus of continuity deteriorates as time 
increases, because of the discontinuity at the boundary (the jump between u and <^ = 0) 
which propagates in the interior. On the contrary, if A (77) < 1, this effect is compensated 
by the fact that the measure /i is not fully integrated, so that we gain an exponential decay 
of the modulus of continuity : some kind of regularization occurs ! (See end Section for 
some examples). 

Let us give two sufficient conditions that explain when both cases may appear (recall 
that n= [j T-^n) : 

Proposition 5.1. // /i(M" \ > 0, then for any i] > 0, 

Kv) < 1 - /"(K" \n) <l. 

On the other hand, if there exists x E such that supp(/x) C T~^fl (with strict inclusion), 
then X{ri) = 1 for 77 > sufficiently small. 

Proof. Let us first notice that obviously, for any x,y e fl, T~^fl fl T~^fl C so that 

fi{T~^Q n Ty^^) < /u(^), which is independent of 77 > and x,y E Q. Thus passing to 
the supremum and using the fact that the total mass of is 1, one gets the first result. 

Now for the second statment: if such an x exists, then for 7/ small enough, any y E Q 
such that |a; — ?/| < rj, will satisfy : supp(/i) C t^^Q fl T~^fl. Hence ij,{t^^Q fl Ty^Q) = 1 
and passing to the supremum for such x, y yields the result. □ 

We refer to the final Section where explicit examples arc given, by let us mention that 
more or less, A > 1 if has a huge support with respect to the shape of fl, whereas A = 1 
if II is rather concentrated (again with respect to the shape of Q). The limit being when 
II — So is so concentrated that nothing moves. 
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Remark 5.1. If if ^ is bounded and continuous, the estimate involves another term 
depending on the integral of ip on A^{ri) = d^Vtj^ \ d^Vt : if 



6(7], t) = max / (p{x + z) d^{z), 



then we can estimate the modulus of continuity using again Gronwall's Lemma, taking 
into account 6{ri,t). The details are left to the reader. 



6. Qualitative Properties of Solutions 
Theorem 6.1. Let us assume that fi > is a measure such that for some rj > 0, 

supp(/i) D {\z\<ri}. 

Then if u{x,t) & X is a solution of problem {V) with if = and uq G C{Q) non-negative 
and compactly supported in Q, we have 

V (x, t) G 1] X (0, oo), u{x, t) > 0. 

Proof. Let us first observe that if we set v{x,t) = e^u{x,t) where u is the extension of u 
defined by fl2.2p . then v satisfies 

dtV = / v{x + z) dfi{z) > , 

so that once v (hence u also) is positive, it stays as such for later times. Let A = 
supp (mo) C Q which is compact and notice first that obviously m > in f2 x [0, oo), using 
the comparison principle. Now, let xq E fl \ A, and assume that dist {xq; A) < rj. Using 
the equation at (xo,0) (remember that the equation holds even on the boundary), one 

gets 

^^t(a;o,0)> / VQ{x + z,0)dfi{z)>0. 

J\z\<rj 

Thus we get that necessarily, for any t > 0, u{x,t) will be positive in = {s G H : 
dist(a;,y4) < r/}. Iterating the same argument at t > arbitrary small, we get that 
u{x,t) > in Akr, for any A; G N. Thus it is clear that u{x,t) > for any x E Cl, for any 
time t > 0. □ 



As a corollary of this result, we have obtained that the boundary data is never taken, 
more precisely : 

Corollary 6.1. In the vanishing viscosity method, a boundary layer occurs and the limit 
solution is strictly positive on dQ for any time t > 0. Thus the boundary data is never 
taken in general. 
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7. Conclusions, Remarks and Examples 

We have proposed a general framework to deal with Dirichlet problems for the non-local 
equation (II .ip . If /i G L^(M"), we have proved existence and uniqueness of a continuous 
solution for any continuous initial and boundary data. Furthermore, under some condition 
on the support of /i near the origin, we have also proved a positivity property, and that 
the boundary values will not be taken in the classical sense. Let us end this note with 
some remarks and extensions : 

1. If /i G L°° , a similar existence and uniqueness result as in Theorem 13.11 may be 
obtained for data UQ,(p in L^, using similar fixed-point arguments. But if /i is a finite 
measure with singular part (like Dirac deltas), then the above construction does not work 
since the integral term may not be continuous. However, we think that similar things 
could maybe be done for measures less concentrated than the delta. 

2. An interesting feature of the equation is that it holds true everywhere on dVt (see 
Proposition 13. ip . This explains why, although the boundary data is not taken in the usual 
sense, there would not be any problem with the viscosity definition of boundary values. 
Indeed, passing to the limit in the viscous equation can be done also in the framework of 
viscosity solutions since the convergence holds in L^^-weak-k. 

3. Theorem 16.11 shows that propagation with infinite speed occurs like for the Heat 
Equation, although the phenomenon is of different nature (due to the non-local term and 
not instantaneous diffusion as it is the case for the heat equation). But even more, u also 
becomes positive at the boundary x G dVL. So, even if both and Uq are continuous, in 
general m(x, t) will not take the boundary data on dVL in the usual sense (this was proved 
for stationary solutions in [5]). And again, we recover the necessity to define a suitable 
notion of boundary for non-local operators, d^Vt here. 

4. Note that some property on the support of /i near the origin is necessary to have 
this propagation property, although our criteria is maybe not optimal. Indeed, if /i is 
too much off-centered, then the operator will not see anything in Vl itself: consider for 
instance =] — 10; +10[, = and let yU be the caracteristic function of the set ] — 1, 0[. 
If u(0) G Co(^^), it will remain compactly supported on the left for any t > 0. 

5. We showed a kind of regularizing effect in the case of bounded domains, provided 
the measure /z satfisfies Xirf) < 1. Let us consider the case f2 =] — 1; 1[ and let us give an 
example of a measure /i satisfying this property : jji{z) = ce"'"^' , where c > is such that 
/i(M) = 1. Then clearly, since Q. =] — 2; 2[, X{r]) < 1 — d/i(2;) < 1. One can also have 
a good estimate of the rate in the exponential decay. Now if ^{z) is the characteristic 
function of the set ] — 1/2; l/2[, we have clearly that for any i] < 1/2, A(?7) = 1. 

6. A different story happens when fi is not integrable near the origin, which is the case 
for the fractional Laplace operator. In this case, the boundary data are taken since the 
measure is so singular that discontinuities at the boundary are not allowed (see [3]). 
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Appendix : The Probabilistic Interpretation 

In this Appendix, we briefly explain the probabiHstic context in which the Cauchy problem 
may be understood. Let (O,^, P) be a measured space and Yq : ^ — ^ 1^ be a random variable 
with law no G L^(]R^) (we assume that J uq = 1). Now let (li)jeN be a sequence of independent, 
identically distributed random variables of law fx, and t i-^ Nt he a, random variable on N with 
Poisson law of intensity t > : 

■p{Nt = k)= e-*- . 

Nt 

Then we consider the stochastic process Xf = Yq + Y^Yj,. The term Yq represents the initial 

1=1 

state of the stochastic process Xt, while the sum represents a random number of jumps Nt, 
each jump being independent of the others, with probability distribution /x. Such a process is 
named a composed Poisson process, which is a special case of the more general Levy processes. 
If wc consider the density u{t) of the probability distribution of Xt, it can be easily shown that 
precisely, 

du 



dt 



I {u{x + z,t) — u{x,t)} djji^z) . 
Jrn 



Now, by definition the characteristic function of Xt is nothing but the Fourier transform 
of u{t), and it is known that for any Levy process, there exists a uniquely determined triplet 
{A, 7, u) such that : 

E(e^^^«) = e*'^(^) with ip{z) = -\z ■ Az + ■ z + [ {e'^^ - 1 - z^lui^iiz)) du{z) . 

In this decomposition, A is a positive-definite matrix, 7 € M" and 1/ is the Levy measure, 
satisfying the following properties : f^^^^-^z^dv , Jj2|>i^^ ^ "^^^ infinitesimal generator of 
the Feller semi-group associated to (Xt) is the following : 

i,j=l i=l 

In the special case of composed Poisson processes, the matrix A is nul, 7 = and the corrector 
term involving the gradient of / vanishes for symmetric measures. 
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